ValueFolder='/sps/isc/rneveu/EtudeIRM/Nature/Results/';
ValueFolderTemp='/sps/isc/rneveu/EtudeIRM/Nature/ResultsJournal/';
ValueFolderTemp2='/sps/isc/rneveu/EtudeIRM/Nature/ResultsEditor/';
load([ValueFolderTemp 'ResultsThirdLevelDistributionOutliers.mat']);
VecDistribution=VecResults;
ValueFolder2='/sps/isc/rneveu/EtudeIRM/Nature/';
load([ValueFolder2 'DonneesNature.mat']);
load([ValueFolder2 'EtudeNature.mat']);
VecEditors1Temp(1,:)=[];
ValueCutoff=1.25;
VecResultsJournalEditor=[];
VecResultsJournalEditorTemp=[];
VecResultsJournalEditorTemp2=[];
VecResultsJournalEditorTemp4=[];
VecResultsJournalEditorTemp3=zeros(2,33);
nEditors=0;
for n2=1:323
    ValueTemp=find(VecDistribution(:,1)==n2);
    if n2~=11&&size(ValueTemp,1)==1&&VecDistribution(ValueTemp(1,1),2)<=ValueCutoff
        ValueTemp=dir([ValueFolderTemp2 'FileAnalysisCharacteristicsEditor' num2str(n2) '.mat']);
        if size(ValueTemp,1)==1
            load([ValueFolderTemp2 'FileAnalysisCharacteristicsEditor' num2str(n2) '.mat']); 
            VecResultsJournalEditor=[VecResultsJournalEditor;VecResultsEditor(1,1:38) n2];
            if strcmp(VecEditors1Temp{n2,13},'Y')==1
                VecResultsJournalEditorTemp2=[VecResultsJournalEditorTemp2;VecResultsEditor(1,1:38) 1];
            else
                VecResultsJournalEditorTemp2=[VecResultsJournalEditorTemp2;VecResultsEditor(1,1:38) 2];
            end
        end
    else
        if n2~=11
            ValueTemp=dir([ValueFolderTemp2 'FileAnalysisCharacteristicsEditor' num2str(n2) '.mat']);
            if size(ValueTemp,1)==1
                load([ValueFolderTemp2 'FileAnalysisCharacteristicsEditor' num2str(n2) '.mat']);
                if strcmp(VecEditors1Temp{n2,13},'Y')==1
                    VecResultsJournalEditorTemp4=[VecResultsJournalEditorTemp4;VecResultsEditor(1,1:38) 1];
                else
                    VecResultsJournalEditorTemp4=[VecResultsJournalEditorTemp4;VecResultsEditor(1,1:38) 2];
                end
            end
        end
    end
    % Robustness analysis including all editors who have data for the duration of research experience
    if n2~=11
        ValueTemp=dir([ValueFolderTemp2 'FileAnalysisCharacteristicsEditor' num2str(n2) '.mat']);
        if size(ValueTemp,1)==1
            load([ValueFolderTemp2 'FileAnalysisCharacteristicsEditor' num2str(n2) '.mat']);
            if strcmp(VecEditors1Temp{n2,13},'Y')==1
                VecResultsJournalEditorTemp=[VecResultsJournalEditorTemp;VecResultsEditor(1,1:38) n2 1];
            else
                VecResultsJournalEditorTemp=[VecResultsJournalEditorTemp;VecResultsEditor(1,1:38) n2 2];
            end
        end
        n3=1;
        TestTemp=0;
        while n3<=33&&TestTemp==0
            if strcmpi(VecJournalTemp{n3,1},VecEditors1Temp{n2,11})==1
                TestTemp=1;
            else
                n3=n3+1;
            end
        end
        if strcmp(VecEditors1Temp{n2,13},'Y')==1
            VecResultsJournalEditorTemp3(1,n3)=VecResultsJournalEditorTemp3(1,n3)+1;
        else
            VecResultsJournalEditorTemp3(2,n3)=VecResultsJournalEditorTemp3(2,n3)+1;
        end
%         load([ValueFolderTemp2 'StartEditor' num2str(n2) '.mat']);
%         if size(ValueMaxDate1,1)==0
%             nEditors=nEditors+1;
%         end
%         if isnan(VecEditors1Temp{n2,3})==1
%             ValueTemp5Temp=VecEditors1Temp{n2,5};
%             if strcmp(ValueTemp5Temp,'NaN')==1||strcmp(ValueTemp5Temp,'No date')==1
%                 if isnan(VecEditors1Temp{n2,22})==1
%                     nEditors=nEditors+1;
%                 end
%             end
%         end 
    end
end

return
ValueTemp22=find(isnan(VecResultsJournalEditor(:,1))==0&isnan(VecResultsJournalEditor(:,27))==0&VecResultsJournalEditor(:,1)<30);
ValueTemp23=find(isnan(VecResultsJournalEditor(:,27))==0&isnan(VecResultsJournalEditor(:,16))==0);
[ValueTemp,ValueTemp2]=corr(VecResultsJournalEditor(ValueTemp22,1)+VecResultsJournalEditor(ValueTemp22,2),VecResultsJournalEditor(ValueTemp22,27),'type','spearman');
[ValueTemp3,ValueTemp4]=corr(VecResultsJournalEditor(ValueTemp22,1)+VecResultsJournalEditor(ValueTemp22,2),VecResultsJournalEditor(ValueTemp22,27),'type','pearson');
[ValueTemp5,ValueTemp6]=corr(VecResultsJournalEditor(ValueTemp23,27),VecResultsJournalEditor(ValueTemp23,16),'type','spearman');
[ValueTemp7,ValueTemp8]=corr(VecResultsJournalEditor(ValueTemp23,27),VecResultsJournalEditor(ValueTemp23,16),'type','pearson');
VecResults=[ValueTemp ValueTemp2 size(ValueTemp22,1);ValueTemp3 ValueTemp4 size(ValueTemp22,1);ValueTemp5 ValueTemp6 size(ValueTemp23,1);ValueTemp7 ValueTemp8 size(ValueTemp23,1)];
ValueTemp23=find(isnan(VecResultsJournalEditorTemp(:,27))==0&isnan(VecResultsJournalEditorTemp(:,16))==0);
[ValueTemp5,ValueTemp6]=corr(VecResultsJournalEditorTemp(ValueTemp23,27),VecResultsJournalEditorTemp(ValueTemp23,16),'type','spearman');
[ValueTemp7,ValueTemp8]=corr(VecResultsJournalEditorTemp(ValueTemp23,27),VecResultsJournalEditorTemp(ValueTemp23,16),'type','pearson');
VecResults=[VecResults;ValueTemp5 ValueTemp6 size(ValueTemp23,1);ValueTemp7 ValueTemp8 size(ValueTemp23,1)];

ValueTemp22=find(isnan(VecResultsJournalEditorTemp2(:,1))==0&isnan(VecResultsJournalEditorTemp2(:,16))==0&VecResultsJournalEditorTemp2(:,1)<30&VecResultsJournalEditorTemp2(:,16)<=nanmedian(VecResultsJournalEditorTemp2(:,16)));
ValueTemp23=find(isnan(VecResultsJournalEditorTemp2(:,1))==0&isnan(VecResultsJournalEditorTemp2(:,16))==0&VecResultsJournalEditorTemp2(:,1)<30&VecResultsJournalEditorTemp2(:,16)>nanmedian(VecResultsJournalEditorTemp2(:,16)));
[ValueTemp,ValueTemp2]=ttest2(VecResultsJournalEditorTemp2(ValueTemp22,1)+VecResultsJournalEditorTemp2(ValueTemp22,2),VecResultsJournalEditorTemp2(ValueTemp23,1)+VecResultsJournalEditorTemp2(ValueTemp23,2));
[ValueTemp3,ValueTemp4]=ranksum(VecResultsJournalEditorTemp2(ValueTemp22,1)+VecResultsJournalEditorTemp2(ValueTemp22,2),VecResultsJournalEditorTemp2(ValueTemp23,1)+VecResultsJournalEditorTemp2(ValueTemp23,2));
VecResults=[VecResults;mean(VecResultsJournalEditorTemp2(ValueTemp22,1)+VecResultsJournalEditorTemp2(ValueTemp22,2)) mean(VecResultsJournalEditorTemp2(ValueTemp23,1)+VecResultsJournalEditorTemp2(ValueTemp23,2)) ValueTemp2;median(VecResultsJournalEditorTemp2(ValueTemp22,1)+VecResultsJournalEditorTemp2(ValueTemp22,2)) median(VecResultsJournalEditorTemp2(ValueTemp23,1)+VecResultsJournalEditorTemp2(ValueTemp23,2)) ValueTemp3];


return

ValueTemp22=find(isnan(VecResultsJournalEditorTemp(:,1))==0&isnan(VecResultsJournalEditorTemp(:,16))==0&VecResultsJournalEditorTemp(:,1)<30&VecResultsJournalEditorTemp(:,40)==2);
ValueTemp23=find(isnan(VecResultsJournalEditorTemp(:,1))==0&isnan(VecResultsJournalEditorTemp(:,16))==0&VecResultsJournalEditorTemp(:,1)<30&VecResultsJournalEditorTemp(:,40)==1);
[ValueResults,ValueResultsTemp,stats]=mnrfit([VecResultsJournalEditorTemp(ValueTemp22,16);VecResultsJournalEditorTemp(ValueTemp23,16)],[ones(size(ValueTemp22,1),1);2*ones(size(ValueTemp23,1),1)]);
[ValueResults2,ValueResultsTemp2,stats2]=mnrfit([VecResultsJournalEditorTemp(ValueTemp22,1)+VecResultsJournalEditorTemp(ValueTemp22,2);VecResultsJournalEditorTemp(ValueTemp23,1)+VecResultsJournalEditorTemp(ValueTemp23,2)],[ones(size(ValueTemp22,1),1);2*ones(size(ValueTemp23,1),1)]);

VecResults3=[];
ValueTemp22=find(isnan(VecResultsJournalEditor(:,1))==0&isnan(VecResultsJournalEditor(:,2))==0);
VecValueTemp=VecResultsJournalEditor(ValueTemp22,2)./VecResultsJournalEditor(ValueTemp22,1);
ValueTemp=find(isnan(VecValueTemp)==1);
VecValueTemp(ValueTemp,:)=ones(size(ValueTemp,1),1);
VecValueTemp=VecValueTemp-1;
[ValueTemp,ValueTemp2]=ttest(VecValueTemp);
VecValueTemp=VecValueTemp+1;
VecResults3=[VecResults3;mean(VecValueTemp*100) sqrt(var(VecValueTemp*100)) ValueTemp2];
ValueTemp22=find(isnan(VecResultsJournalEditor(:,1))==0&isnan(VecResultsJournalEditor(:,24))==0&isnan(VecResultsJournalEditor(:,25))==0);
VecValueTemp=(VecResultsJournalEditor(ValueTemp22,24)+VecResultsJournalEditor(ValueTemp22,25))./VecResultsJournalEditor(ValueTemp22,3);
ValueTemp=find(isnan(VecValueTemp)==1);
VecValueTemp(ValueTemp,:)=ones(size(ValueTemp,1),1);
VecValueTemp=VecValueTemp-1;
[ValueTemp,ValueTemp2]=ttest(VecValueTemp);
VecValueTemp=VecValueTemp+1;
[ValueTemp3,ValueTemp4]=ttest(VecResultsJournalEditor(ValueTemp22,24)+VecResultsJournalEditor(ValueTemp22,25),VecResultsJournalEditor(ValueTemp22,3));
VecResults3=[VecResults3;mean(VecValueTemp*100) sqrt(var(VecValueTemp*100)) ValueTemp2;mean(VecResultsJournalEditor(ValueTemp22,24)+VecResultsJournalEditor(ValueTemp22,25)) sqrt(var(VecResultsJournalEditor(ValueTemp22,24)+VecResultsJournalEditor(ValueTemp22,25))) ValueTemp4];
ValueTemp22=find(isnan(VecResultsJournalEditor(:,1))==0&isnan(VecResultsJournalEditor(:,2))==0&VecResultsJournalEditor(:,1)<30)
VecValueTemp=VecResultsJournalEditor(ValueTemp22,2)./VecResultsJournalEditor(ValueTemp22,1);
ValueTemp=find(isnan(VecValueTemp)==1);
VecValueTemp(ValueTemp,:)=ones(size(ValueTemp,1),1);
VecValueTemp=VecValueTemp-1;
[ValueTemp,ValueTemp2]=ttest(VecResultsJournalEditor(ValueTemp22,1),VecResultsJournalEditor(ValueTemp22,2))
VecResults3=[VecResults3;mean(VecResultsJournalEditor(ValueTemp22,1)) sqrt(var(VecResultsJournalEditor(ValueTemp22,1))) ValueTemp2;mean(VecResultsJournalEditor(ValueTemp22,2)) sqrt(var(VecResultsJournalEditor(ValueTemp22,2))) NaN];
ValueTemp22=find(isnan(VecResultsJournalEditorTemp2(:,1))==0&isnan(VecResultsJournalEditorTemp2(:,11))==0&VecResultsJournalEditorTemp2(:,1)<30&VecResultsJournalEditorTemp2(:,39)==2);
VecResults5=[mean(VecResultsJournalEditorTemp2(ValueTemp22,11)) sqrt(var(VecResultsJournalEditorTemp2(ValueTemp22,11)))];
ValueTemp22=find(isnan(VecResultsJournalEditorTemp2(:,1))==0&isnan(VecResultsJournalEditorTemp2(:,16))==0&VecResultsJournalEditorTemp2(:,1)<30&VecResultsJournalEditorTemp2(:,39)==2);
ValueTemp23=find(isnan(VecResultsJournalEditorTemp2(:,1))==0&isnan(VecResultsJournalEditorTemp2(:,16))==0&VecResultsJournalEditorTemp2(:,1)<30&VecResultsJournalEditorTemp2(:,39)==1);
[ValueTemp,ValueTemp2]=ttest2(VecResultsJournalEditorTemp2(ValueTemp22,16),VecResultsJournalEditorTemp2(ValueTemp23,16));
VecResults3=[VecResults3;mean(VecResultsJournalEditorTemp2(ValueTemp22,16)) sqrt(var(VecResultsJournalEditorTemp2(ValueTemp22,16))) ValueTemp2];
VecResults3=[VecResults3;mean(VecResultsJournalEditorTemp2(ValueTemp23,16)) sqrt(var(VecResultsJournalEditorTemp2(ValueTemp23,16))) NaN];

% Comparison of the duration of the Phd and the number of original articles published during the PhD between editors with only a Phd and the others
VecResults2=[];
ValueTemp22=find(isnan(VecResultsJournalEditorTemp(:,12))==0);
VecValueTemp=VecResultsJournalEditorTemp(ValueTemp22,:);
ValueTemp=find(VecValueTemp(:,end)==1);
ValueTemp2=find(VecValueTemp(:,end)==2);
[ValueTemp3,ValueTemp4]=ttest2(VecValueTemp(ValueTemp,12),VecValueTemp(ValueTemp2,12));
VecResults2=[VecResults2;mean(VecValueTemp(ValueTemp,12)) sqrt(var(VecValueTemp(ValueTemp,12))) mean(VecValueTemp(ValueTemp2,12)) sqrt(var(VecValueTemp(ValueTemp2,12))) ValueTemp4];
ValueTemp22=find(isnan(VecResultsJournalEditorTemp2(:,11))==0&VecResultsJournalEditorTemp2(:,1)<30&VecResultsJournalEditorTemp2(:,end)==2);
ValueTemp23=find(isnan(VecResultsJournalEditorTemp4(:,11))==0&VecResultsJournalEditorTemp4(:,end)==2);
ValueTemp24=find(isnan(VecResultsJournalEditorTemp2(:,11))==0&VecResultsJournalEditorTemp2(:,1)<30&VecResultsJournalEditorTemp2(:,end)==1);
ValueTemp25=find(isnan(VecResultsJournalEditorTemp4(:,11))==0&VecResultsJournalEditorTemp4(:,end)==1);
[ValueTemp3,ValueTemp4]=ttest2(VecResultsJournalEditorTemp2(ValueTemp22,11),VecResultsJournalEditorTemp4(ValueTemp23,11));
VecResults2=[VecResults2;mean(VecResultsJournalEditorTemp2(ValueTemp22,11)) sqrt(var(VecResultsJournalEditorTemp2(ValueTemp22,11))) mean(VecResultsJournalEditorTemp4(ValueTemp23,11)) sqrt(var(VecResultsJournalEditorTemp4(ValueTemp23,11))) ValueTemp4];
[ValueTemp3,ValueTemp4,stats]=fishertest([size(ValueTemp22,1) size(ValueTemp23,1);size(ValueTemp24,1) size(ValueTemp25,1)])
return
ValueTemp22=find(isnan(VecResultsJournalEditorTemp2(:,11))==0);
VecValueTemp=VecResultsJournalEditorTemp2(ValueTemp22,11);
VecResults2=[VecResults2;mean(VecValueTemp) sqrt(var(VecValueTemp)) NaN NaN NaN];
ValueTemp22=find(isnan(VecResultsJournalEditorTemp2(:,13))==0);
VecValueTemp=VecResultsJournalEditorTemp2(ValueTemp22,:);
ValueTemp=find(VecValueTemp(:,end)==1);
ValueTemp2=find(VecValueTemp(:,end)==2);
[ValueTemp3,ValueTemp4]=ranksum(VecValueTemp(ValueTemp,13),VecValueTemp(ValueTemp2,13));
VecResults2=[VecResults2;mean(VecValueTemp(ValueTemp,13)) sqrt(var(VecValueTemp(ValueTemp,13))) mean(VecValueTemp(ValueTemp2,13)) sqrt(var(VecValueTemp(ValueTemp2,13))) ValueTemp4];

ValueResult=VecResultsJournalEditorTemp3(1,32)/(VecResultsJournalEditorTemp3(1,32)+VecResultsJournalEditorTemp3(2,32))*100;
ValueTemp=1:31;
ValueTemp=[ValueTemp 33];
[ValueTemp2,ValueTemp3]=ttest(VecResultsJournalEditorTemp3(1,ValueTemp)./(VecResultsJournalEditorTemp3(1,ValueTemp)+VecResultsJournalEditorTemp3(2,ValueTemp))*100-ValueTemp);
VecResults4=[ValueResult mean(VecResultsJournalEditorTemp3(1,ValueTemp)./(VecResultsJournalEditorTemp3(1,ValueTemp)+VecResultsJournalEditorTemp3(2,ValueTemp))*100) sqrt(var(VecResultsJournalEditorTemp3(1,ValueTemp)./(VecResultsJournalEditorTemp3(1,ValueTemp)+VecResultsJournalEditorTemp3(2,ValueTemp))*100)) ValueTemp3];
% save([ValueFolderTemp 'AnalyseCorrelationEditors.mat'],'VecResultsJournalEditor');
